Optimized single-qubit gates for Josephson phase qubits 
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In a Josephson phase qubit the coherent manipulations of the computational states are achieved 
by modulating an applied ac current, typically in the microwave range. In this work we show that it 
is possible to find optimal modulations of the bias current to achieve high-fidelity gates. We apply 
quantum optimal control theory to determine the form of the pulses and study in details the case 
of a NOT-gate. To test the efficiency of the optimized pulses in an experimental setup, we also 
address the effect of possible imperfections in the pulses shapes, the role of off-resonance elements 
in the Hamiltonian, and the effect of capacitive interaction with a second qubit. 

PACS numbers: 85.25.Cp, 03.67.Lx, 02.30.Yy 



I. INTRODUCTION 



Over the past decades, together with the development 
of the theory of quantum information— there has been 
an increasing effort to find those physical systems where 
quantum information processing could be implemented. 
Among the many different proposals, devices based on 
superconducting Josephson junctions are promising can- 
didates in the solid state realm (see the reviews in Ref. 
0111111). Josephson qubits can be categorized into three 
main classes: Charge, phase and flux qubits, depending 
on which dynamical variable is most well defined and con- 
sequently which basis states are used as computational 
states |0) and 

Phase qubits^iii^, subject of the present investigation, 
in their simplest configuration can be realized with a sin- 
gle current biased Josephson junction. For bias lower 
that the critical current, the two lowest eigenstates of 
the system form the computational space. The appli- 
cation of a current pulse, with frequency which is in 
resonance with the transition frequency of the two log- 
ical states, typically in the microwave range, allows to 
perform all desired single bit operations. Recent experi- 
ments^iiS have realized both single-bit and two-bit gates 
in capacitive coupled phase qubits. In the experiments 
conducted so far, motivated by similar approach in NMR, 
the amplitude of the microwave current used to per- 
form the qubit manipulation has a Gaussian shape^i^i. 
The importance of achieving fast quantum gates with 
high fidelity rises the question whether there are mod- 
ulations, other than Gaussian, which lead to higher fi- 
delities. Indeed, Gaussian pulses act better when their 
duration time is longer—, therefore the search for modu- 
lations which result in high fidelity gates, even when the 
duration time of the pulse is short, seems to be neces- 
sary. Some theoretical work has already been done in 
this direction to examine the efficiency of different mod- 
ulationsii. In the present paper we follow a different 
approach as compared to Ref. |ll| and show that by em- 
ploying the quantum optimal control theorjsi^ii^iii^ii^, we 
can further improve the (theoretical) bounds on the error 



of gate operations. 

Quantum optimal control has been already applied to 
optimize quantum manipulation of Josephson nanocir- 
cuits in the charge limifei^iiiii^. Here we want to test 
this method in the opposite regime of phase qubilii^ and 
see whether it is possible to find optimal modulations of 
microwave pulses, with different duration times, which 
give very good fidelity for single bit operations. 

The paper is organized as follows: in Sec. |TT] we will 
describe the model for the phase qubit and the Hamil- 
tonian used in the rest of the paper. In Sec. IIIII we in- 
troduce the NOT quantum gate which we have chosen 
to optimize. Then a brief introduction to the quantum 
optimal control algorithm which is used for this work will 
be given in Sec. IIVI The numerical results for a phase 
qubit will be presented in Sec. El The achieved accu- 
racy for desired operation, discussed in Sec. IV Al is fur- 
ther tested against possible imperfections in the pulses 
shape (Sec. IV B[) . presence of off- resonance elements in 
the Hamiltonian (Sec. IV Cp and possible presence of the 
inter-qubit capacitive interaction in multi-qubit systems 
(Sec lVDj) . The specific question of the leakage out of the 
Hilbert space is addressed in Sec. IVII where we provide 
numerical results obtained for a junction with five levels 
inside its potential. A summary of the results obtained 
and possible perspectives of this work will be presented 
in the concluding remarks in Sec. IVIll 



II. SINGLE-JUNCTION PHASE QUBIT 

A phase qubit can be realized by a fiux-biased rf 
SQUII>2S, a low inductance dc SQUID^ or a large in- 
ductance dc SQUID^. In its simplest design a phase 
qubit consists of a single Josephson junction (Fig. [H^a)) 
with critical current Iq and a biasing dc current I^c- The 
Hamiltonian has the form 



(1) 



where Ec ~ (2e)^/2C and Ej = /o$o/27r are, respec- 
tively, the charging energy and the Josephson energy of 
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(a) (b) 

FIG. 1; (a) Schematic drawing of a single-junction phase 
qubit with capacitance C, resistance R and critical current 
Jo which is biased by a dc current I^c- A microwave pulse 
Itiw ~ I{t) cos(t<;t + (fi), with frequency u — loqi, is applied to 
make transitions between the two lowest energy levels of the 
system |0) and (b) By neglecting the resistive branch, 
the potential energy of the system U, as a function of the 
Josephson phase across the junction S, has the form of a tilted 
washboard potential. This potential is defined by the height 
of the well AU and the frequency of the classical oscillations 
in the bottom of the well. 



the junction with capacitance C, $o — h/2e being the 
quantum of flux, and S represents the Josephson phase 
across the junction. The regime in which the supercon- 
ducting phase S is the appropriate quantum variable is 
reached when Ej ^ Ec- The potential energy of the sys- 
tem, as a function of 6, has the form of tilted washboard 
with quantized energy levels inside each well (Fig.[TJb)). 
When Idc ^ ^0 there are few levels inside each well and 
the two lowest states |0) and with energies Eq and 
El and transition frequency wqi = (-^'i ~ E^j/hK, 5 GHz, 
can be used as computational states. The transition 
between the two lowest states is made by use of a mi- 
crowave current /^^ of frequency tn which is in resonance 
with transition frequency wqi • Transition to higher states 
(|2), |3) ...), which are out of qubit manifold, are off- 
resonance due to the anharmonicity of the potential well. 

To be effectively used as a two-level quantum system, 
the junction is biased with a dc current slightly smaller 
than the critical current Idc ^ lo- In this regime the 
potential energy of the system can be approximated by 
a cubic potential and the Hamiltonian ([T|) becomes: 

-^(^-f)^ (2) 

The application of a microwave current I^^j = 
I{t) cos{ujt + ip) is taken into account by adding the linear 
term iJ^u, = to the Hamiltonian ([1]). Since the 

eigenstates of the junction biased with a dc current are 
used as computational states, it is appropriate to write 
the full Hamiltonian in the basis of the eigenstates |n) of 
the system with dc bias current. To examine the effect 
of microwave current, one needs to know the elements of 



the superconducting phase S in this basis. 

Moving to the rotating frame, in which the fast oscilla- 
tions due to cos{Lut + (f) do not appear, the Hamiltonian 
H in the rotating frame is related to the Hamiltonian in 
laboratory frame H via 

H = VHV'' -ihV^V'', (3) 

whereas the state of the system in the rotating frame is 
l^/i) = V\i}). By introducing g{t) = I {t) ^/h/2(Mn and 
Amn = ^Vf^oiiiTT-lSln) , and considering only the first 
three levels in the well, the Hamiltonian of the phase 
qubit in the rotating frame takes the following form 

/ 5(t)Aoie^'^ \ 

H « .g(i)Aioe-*'^ g{t)Aue'^ . (4) 

\ g{t)A2ie-"^ -hSuj J 

Here we have set £'o = 0, to = loqi, Soj = luqi — 0^12 and 
we have assumed that off-resonance terms have negligible 
effect. As we shall see, by a proper choice of if and mi- 
crowave current modulation g{t) it is possible to perform 
single-bit operations on the computational states |0) and 

III. NOT-GATE 

As one can see from the 2x2 top-left block of the 
Hamiltonian ([5]), the initial phase of the microwave pulse 
if defines the axis of rotation, in xy-plane of the Bloch 
sphere, for a given state, while the pulse amplitude and 
duration time define the angle of rotation. For example, 
by setting (p = {(p = tt/2) such block is proportional 
to the Pauli matrix ax (cry), i. e. a rotation around the 
x(y)-axis. In a recent experiment^ a tt rotation around 
X has been implemented as a part of a sequence of oper- 
ations to create entanglement between two phase qubits. 
This motivates us to set ip = and focus this work on 
the single-qubit NOT-gate operation consisting of a tt 
rotation around the x-axis. 

In the typical experiment a shaped pulse with the fol- 
lowing Gaussian modulationii 

g{t) ^ - e (5) 
'-a 

is used to induce flips between states |0) and |1) and 
vice versa. Here a, tg and T = 2atg are, respectively, 
the amplitude, characteristic width and total width of 
the pulse, a being the cut-off of the pulse in time. The 
actual result of the operation can be quantified by the 
fidelity \{ip{T)\ipfin)\'^ , where \ipfin) is the desired final 
state and \tp{T)) is the state achieved at the end of time 
evolution starting from initial state \i}!{t = 0)) = \ipini)- 

For a TT rotation and with a typical cut-off value 
(3 < a < 5) the amplitude a w \pnj2 yields a pretty 
high fidelity of rotation. More precisely. Fig. [5] shows the 
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FIG. 2: (Color on line) Error E for a NOT-gate operation on 
an arbitrary superposition of states |0) and |1) after applying 
a Gaussian pulse with amplitude a = 1.25, cut-off a = 3 as a 
function of the duration time T. 




FIG. 3: (Color on line) The leakage outside the qubit mani- 
fold for a NOT-gate operation on an arbitrary superposition 
of states |0) and |1) after applying a Gaussian pulse with am- 
plitude a — 1.25, cut-off a = 3 as a function of the duration 
time T. 



error £ = 1 — |(V'(r)|'0/m) P for a NOT-gate operation 
on an arbitrary superposition jV'irii) = b\0) + c|l), which 
would result in the state \ipfin) = + c|0), using a 
Gaussian pulse with cut-ofF a = 3, amplitude a = 1.25 
and duration time T. The leakage outside the qubit man- 
ifold, defined as \{i:{T)\2)\'^ , is also shown in Fig. El It 
is worthwhile noting that although the leakage, for long 
enough pulses, can be of the order of 10~^, the error in 
the NOT-gate operation is always higher than 10~^ (see 
Ref.lHl). 



IV. QUANTUM OPTIMAL CONTROL 

As we mentioned in the Introduction, in this work we 
use quantum optimal control theory in order to obtain 
microwave current modulations which give rise to a high- 
fidelity NOT-gate operation for a phase qubit. In this 
section we briefiy review the optimal control algorithm 



which we have employed to obtain optimized modula- 
tions. 

In general, quantum optimal control algorithmsi^J^ii^ 
are designed to lead a quantum system with state \ip{t)) 
from an initial state |V'(0)) = IV'ini) to a target final state 
\ipfin) at time T by minimizing a cost functional which 
is a measure of inaccuracy of reaching the desired final 
state. If \ip{T)) denotes the state achieved at time T, one 
can consider two different cost functionals: 

. ei = l-|(V'(T)|^/,„)P 
By minimizing this cost functional, although the 
population of the desired state I'lpfin) will be max- 
imized, the overall phase of this state is not forced 
to be preserved. 

Minimization of this second cost functional, in ad- 
dition to maximizing the population of the desired 
state, preserves its overall phase. 

In optimal control theory the minimization of the cost 
functional is done by updating the Hamiltonian of the 
system, via some control parameters, in an iterative pro- 
cedure until the desired value of the cost functional is 
reached. Any specific algorithm which is guaranteed to 
give improvement at each iteration^^ is called immediate 
feedback control and can be briefly described as follows: 
Assume that the Hamiltonian of the system depends on 
a set of parameters {uj{t)} which are controllable. By 
using a proper initial guess {u^^\t)} for control param- 
eters, first the state of the system \ip{t)) is evolved in 
time with the initial condition \ip{0)) = IV'ini) giving rise 
to lipiT)) after time T. At this point the iterative algo- 
rithm starts, aiming at decreasing the cost functional by 
adding a correction to control parameters in each step. 
In the nth step of this iterative algorithm 

• An auxiliary state |x(0) is evolved backward in 
time starting from \x(T)) reaching |x(0)). 

In the case of minimizing ei, IxC^)) = 
\'<Pfin){'4'fin\tpiT)) and for minimizing 62, \x{T)) = 

2mT)) - l^f^r.)). 

• The states |x(0)) and |'0(O)) are evolved for- 
ward in time, respectively, with control parameters 
{u^"^(t)} and {u^"+'\t)}. Here, 



(ri+l)/,N {") u\ I c 



{x{t)\ 



dH 

duj {t) 



\m) 



(6) 



are updated control parameters. X(t) is a weight 
function used to fix initial and final conditions on 
the control parameters in order to avoid major 
changes at the beginning and end of time evolution 
and is an important parameter for the convergence 
of the algorithm. 
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These two steps are repeated until the desired value of 
ei or 62 is obtained. 

In order to implement the optimization procedure to 
a NOT-gate for any arbitrary superposition of compu- 
tational states, one must be able to flip |0) and |1) at 
the same time (i. e. with same pulse) making sure that 
the phase relation between them is preserved. This is 
guaranteed by using the following definition of fidelity 



T = 



(7/.o(T)|l) + (V.i(T)|0) 



(7) 



where |'(/'o(T')) and |-!/'i(r)) are final states achieved at 
time T after applying the same pulse on initial states 
|0) and |1). The minimization of the cost functional ei, 
for flipping at the same time the states |0) and does 
not necessarily lead to maximization of the fldelity ([7]) 
due to possible changes in the phase relation between 
them. However if ei is minimized, the maximal fldelity 
is also guaranteed. Therefore in order to obtain a high- 
fidelity NOT-gate it seems more natural to minimize 62 
instead of e\. However, in the following we will show 
that, although in the ideal case optimized pulses obtained 
from minimizing 62 result in much higher fidelity, when 
more realistic cases are considered optimized pulses from 
minimizing 61 lead to higher fidelities, specially for very 
short pulses. In this work we often use the error £ ^ \—T 
instead of fidelity. 



V. NUMERICAL RESULTS 

In this section we present the numerical results to show 
that the quantum optimal control theory allows to op- 
timize the modulation of microwave pulses in order to 
implement a high-fidelity NOT-gate. The optimization 
is done in the rotating frame and the Hamiltonian ([4]) is 
used for time evolution while Aj^ are calculated by means 
of perturbation theory. 



A. Optimal NOT-gate 

By employing the quantum optimal control algorithm 
described in the Sec. IIVI and using the modulation of 
the microwave pulse g(i) as the control parameter, we 
start from Gaussian pulses (O of given duration time T 
as the initial guess and optimize the NOT-gate opera- 
tion. We will show the results obtained from minimizing 
both 61 and 62 and refer to corresponding errors by £\ 
and El and corresponding optimized pulses by g\ and gi . 
The optimization has been stopped when either the cost 
functionals reached the value 10~^^ or 5000 iterations are 
done. 

Figure [4] shows the error f , as a function of duration 
time of the pulse T, for the Gaussian pulses used as initial 
guess (circles) and for the optimized pulses (unfilled tri- 
angles and squares). For most of points, the convergence 




FIG. 4: (Color on line) Error for a NOT-gate operation ap- 
plied to any arbitrary superposition of states |0) and |1) made 
by pulses with Gaussian modulation (circles) and optimized 
modulation obtained from minimizing ei (unfilled triangles) 
and £2 (squares) in a three-level system. Gaussian pulses 
have amplitude a — 1.25 and cut-off in time a = 3. Opti- 
mized pulses are obtained after at most 5000 iterations and 
using Gaussian pulses as initial guess. Filled triangles are ob- 
tained by applying a O.OItt phase shift after optimized pulses 
obtained from minimizing ei. 



is reached in much less than 5000 iterations. However 
for pulses with T < 2 5000 iterations has been com- 



Si^' , . . . 

plcted. As we expected, minimizing 62 results in high- 
fidchty NOT-gate with £ « 10" for aU T > 2 |^ w 4 
ns, while for very short pulses it seems that, with same 
number of iterations, minimizing 61 leads to better re- 
sults. 

In order to understand the reason for the oscillating 
behavior of £1 as a function of T, we plot the average 
value of 61 for |G) |1) and |1) |0) transitions at the 
end of optimization (top panel of Fig. [5]) which shows 
that the final value of 61 for both of these transitions is 
of the order of 10~^^. As we explained in Sec. IIVI 61 
is insensitive to the phase of the final state and it turns 
out that while a given optimized pulse applied to initial 
state |0) leads to the final state 6**''|1), then the same 
pulse might transform the initial state |1) into 6**i|0), i. 
e. there is a phase difference between the two final states 
Odiff = 01 — 00- The bottom panel of Fig. [5] shows this 
phase difference for optimized pulses with given duration 
time T which increases the error of NOT-gate £1 to what 
has been shown in Fig. 21 

Although this phase difference causes a major increase 
in the error while working with superpositions, the error 
£1 is at least one order of magnitude smaller than those 
from Gaussian pulses (Fig. [3]) . Moreover the final phase 
difference between |0) and |1) can be compensated by a 
following phase shift gate. In Fig.|3]results after applying 
a O.OItt (which is approximately the average of 0diff in 
time) phase shift are also shown (filled triangles) which 
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FIG. 5: (Color on line) Top panel: the averaged value of ei 
for |0) <-> |1) transitions after applying pulses with Gaussian 
modulation (circles) and optimized modulation (triangles). 
Optimized pulses are obtained by minimizing ei. Vertical 
a:xis is in logarithmic scale and T is the total width of the 
pulse. Bottom panel: the phase difference between final |0) 
and |1) states (after applying the optimized pulses) in units 
of TT. This final phase difference increases the error of NOT- 
gate £i to what has been shown in Fig. |4] In principle, a 
proper phase shift gate can compensate this phase difference 
and decrease the error to fO~^^. 

demonstrate a significant decrease of Ei. 

B. Imperfections in the pulse shapes 

In this section we study the Fourier transform of the 
optimized pulses, in order to see how practically they are 
realizable in the laboratory, and to examine the effect of 
high-frequency components. Two examples of the final 
optimized pulses (dashed lines) are shown in Fig.[6l both 
with duration time T = 2 |^ w 4 ns. Optimized gi (top 
panel) and (72 (bottom panel) are the results of minimiz- 
ing, respectively, ei and 62 which for T = 2 |^ both arc 
of the order of 10~^^. The corresponding Gaussian pulse 
is also shown in both panels, gi is guaranteed to decrease 
the error of NOT-gatc two orders of magnitude with re- 
spect to the Gaussian pulse while 52 would reduce the 
error up to ten orders of magnitude. 

Fig.[7]shows the Fourier transform of the two optimized 
pulses shown in Fig. [H] To filter out the high-frequency 
components of the optimized pulses we set a cutoff fre- 
quency oJcut for Fourier components and apply the trun- 
cated pulses again and obtain the error. Fig.[5]shows the 
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FIG. 6: (Color on line) Examples of final optimized modula- 
tion of pulses (dashed lines) and the corresponding pulse with 
Gaussian modulation (solid lines) used as initial guess in opti- 
mization process with duration time T = 2 |^ obtained from 
minimizing ei (top panel) and 62 (bottom panel). 



error for a NOT-gate for pulses with different duration 
times as functions of ujcut- t^oi/27r is approximately 5 
GHz and 5uj is typically 10% of woi- In our calculation 
Slo = 0.1 woi which means that 5lo /2ti 500 MHz. 

In the case of £1 , top panel of figure [8] makes clear 
that all important harmonics have frequencies smaller 
than ^5lo. Note that the number of harmonics in- 
cluded within the cutoff is equal to Tujcutl'^''^ so that, 
for T = N{2tt/Slli), such number is equal to N times the 
ratio LOcutl^^- As a result it seems that, for all values 
of T considered, about 20 harmonics should be sufficient 
to reach the smallest value of £\. £1, though, seems to 
be more sensitive to high-frequency components but still 
about four orders of magnitude smaller than £\ under 
the cutoff LOcut = lOJw. 



C. Effect of off-resonance terms 

As we mentioned before, we have assumed that off- 
resonance elements of the Hamiltonian ([3]) in the rotat- 
ing frame are negligible and we have used Hamiltonian 
dH) for calculating the evolution. In this section we check 
this assumption by addressing the effect of off-resonance 
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FIG. 7: (Color online) Fourier transform g{uj) of two opti- 
mized pulses shown in Fig. [6] (dashed lines). gi{ui) minimizes 
ei and (72(1^) minimizes 62. (5tJ is chosen to be ten percent of 
ujoi and [.tJoi/27r is approximately 5 GHz. 
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elements by evolving the complete Hamiltonian ^ us- 
ing the optimized pulses obtained using Hamiltonian ([4]). 
Top panel of Fig. [9] shows the error for a NOT-gate oper- 
ation implemented by Gaussian (circles) and optimized 
pulses from minimizing ei (triangles) and 62 (squares). 
For T > 2 1^ the optimized pulses yield a much higher 
error, with respect to the case when off-resonance terms 
are neglected, still showing an improvement of two or- 
ders of magnitude if compared to Gaussian pulses. Bot- 
tom panel of Fig. [9] shows the absolute value of the er- 
ror difference 6£ obtained by subtracting the error with- 
out off-resonance term from the error with off-resonance 
terms. These figures make clear that while for Gaus- 
sian pulses off-resonance terms can be neglected, for op- 
timized pulses, specially those obtained from minimizing 
62, they are very important. Note that, contrary to the 
ideal case where £2 was about eight orders of magnitude 
smaller than £1, under the effect of off-resonance terms, 
£2 seems to be larger than £1 specially for very short 
pulses with T < 2 This means that the assump- 
tion of ignoring these terms is more accurate when ei is 
minimized. The simpler shape of the optimized pulses 
obtained from minimization of ei could be a reason for 
that. 



FIG. 8: (Color on line) Error £ for a NOT-gate with opti- 
mized pulses obtained by minimizing ei (top panel) and 62 
(bottom panel) as a function of the cutoff frequency ujcut- In- 
teger values of Tujcut/'i'T correspond to the number of Fourier 
components included. 5Lu/2n is approximately 500 MHz. 



D. Effect of capacitive interaction 

So far we have considered a single qubit with three en- 
ergy levels and obtained the modulation of the microwave 
pulses in order to optimize the NOT-gate operation for 
the two lowest energy states |0) and It is now in- 
teresting to consider the setupJ^ containing two qubits 
interacting via a capacitor. The question that we want 
to address is what happens if these optimized pulses are 
applied on the first qubit while the interaction with the 
second qubit is present. 

The interaction Hamiltonian of a circuit with two iden- 
tical phase qubits has the following form: 

where Si and 62 are Josephson phases across the junc- 
tion 1 and 2 and Cx is the capacitance of the interaction 
capacitor. Note that the term with second derivative 
in Eq. ([5]) can be included in the Hamiltonians of the 
uncoupled qubits ([1]) by replacing the charging energy 
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FIG. 9: (Color on line) Top panel: the error £ for a NOT- 
gate made by applying microwave pulses with Gaussian mod- 
ulation (circles) and optimized modulation (triangles and 
squares) when ofF-resonance terms are kept. Note that op- 
timized pulses are obtained by excluding off-resonance ele- 
ments. Bottom panel: the absolute value of error difference 
S£ obtained by subtracting the curves in the top panel from 
those in Fig. [J] T is the total time width of the pulses. 

Ec with an eflFective one £^Ce// — (2e)^/(2Ce//), where 
Ceff = C^/Cs and Ce = C + C^. The Hamiltonian 
can again be written in the basis of the cigenstates of 
the uncoupled qubits, and the strength of the interaction 
Hamiltonian reduces to (Cx/Cs){huJoi). We move to the 
rotating frame described by the unitary operator 

/I \ /I \ 

V = [O e*'^* (8) e'"* (9) 
\0 e^'^^y \0 e2*'^7 

and neglect the off-resonance elements of the result- 
ing Hamiltonian. By applying microwave pulse on the 
first qubit, our aim is to perform a NOT-gate opera- 
tion on such qubit (namely. (TxI ® l2)- Since is 
typically of few fF and C is of the order of pF— we 
find Cx/Cs ~ 2.3 x 10^'^ which leads to an interaction 
strength {Cx/C^)uJoi ~ lOMHz. Figure [TU] shows the 
error as function of time width of the pulse T for both 
Gaussian (circles) and optimized pulses ( triangles and 
squares). Although the optimized pulses are obtained 
for a single qubit system, they still result in smaller error 
at least for short pulses. These results show the impor- 
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FIG. 10: (Color on line) Error £ for a NOT-gate on the first 
qubit implemented by applying Gaussian (circles) and opti- 
mized (triangles and squares) pulses on the first qubit of a 
two-qubit system in presence of capacitive interaction. Gaus- 
sian pulses have amplitude a = 1.25 and cutoff in time a — 3 
and are used as initial guess in the optimization procedure. 
Optimized pulses are obtained from the single-qubit setup. 
The strength of the interaction is Cx/Cs = 2.3 x 10~^. Bot- 
tom panel: the absolute value of error difference 5£ obtained 
by subtracting the curves in the top panel from those in Fig.|4l 



tancc of the presence of the capacitive interaction even 
though the strength of the interaction is small. As it is 
clear from Fig. [TO] for pulses longer than, approximately, 
8 ns (T = 4 1^) the interaction becomes more effective 
and the error for Gaussian and optimized pulses are very 
close. Moreover, longer pulses lead to higher value of er- 
ror contrary to what happens in the case of a single qubit. 
In this case also for very short pulses (T < 1.75 |^) £i 
is smaller than £2 and for T > 2 |^ they are of the same 
order, although in the ideal case £1 was eight orders of 
magnitude larger than £2- 

Bottom panel of Fig. [TO] shows the absolute value of 
error difference S£ which is obtained by subtracting the 
error of ideal case from the error in presence of interac- 
tion. S£ is almost the same in all three cases. 



VI. LEAKAGE 

As explained in section [Hi the two lowest energy levels 
of a current-biased Josephson junction can be used as |0) 
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FIG. 11: (Color on line) The error £ of NOT-gate made 
by applying Gaussian (circles) and optimized (triangles and 
squares) pulses as function of duration time of the pulse T. 
Gaussian pulses have amplitude a = 1.25 and cutoff in time 
a = 3 and are used as initial guess in optimization. Optimized 
pulses are obtained after at most 15000 iterations. Physical 
system contains five energy states and Suj is assumed to be 
0.05 tJoi. 
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FIG. 12: (Color on line) Top panel: the averaged value of ei 
for |0) <^ |1) transitions after applying the Gaussian pulses 
(circles) and optimized pulses (triangles) obtained from mini- 
mizing ei. Vertical axis is in logarithmic scale. Bottom panel: 
the phase difference between final |0) and |1) states in units 

of TT. 



and |1) states of a phase qubit. Although it would be 
desirable to have only two levels inside the potential well 
of the Fig. [1] this is not the case in experimental setups. 
So far we have included the leakage by considering only 
an additional third level and showed that it is possible 
to optimize the pulses in order to gain high fidelity for a 
NOT-gate for a single qubit. In typical experiments the 
number of energy levels inside the well varies between 
three and five. In order to have a more complete under- 
standing of the leakage, in this section we show some re- 
sults obtained for a five-level system. Since adding more 
levels to the system decreases the inhomogencity of the 
level-spacing we choose 5uj — O.OSwqi. 

Figure [11] shows the error £ for a NOT-gate imple- 
mented by optimized and Gaussian pulses, which are 
used as initial guess, for different duration times T . Sim- 
ilar to the case of three-level system, with same number 
of iterations, minimization of ei leads to better results 
for short pulses while for longer duration times of pulses 
minimizing 62 results in error of NOT-gate £2 ~ 10~^^. 
In the case of minimizing ei, at least one order of mag- 
nitude improvement is achieved for long pulses, although 
the improvement obtained for pulses with shorter time 
width are the best. By looking at the average value of ei 
for transitions between the states |0) and |1) (top panel 
of Fig. [T^ and the final phase difference between them 
(bottom panel of Fig. [72]) . one realizes that, as it was ob- 
served in three-level system, considerable amount of £1 
is due to the final phase difference Odiff- For instance 
the pulse with T = 2 |^ results in a phase difference ap- 
proximately equal to zero and therefore £1 for this pulse 
is of the order of 10~^. A proper phase shift applied af- 
ter the NOT-gate operation will compensate the phase 
difference between the final |0) and |1) states and conse- 



quently attaining a very high fidelity. 

Two examples of pulses with T = 2 {2-k /5lj) are shown 
in figure [ini gi is obtained from minimizing ei and gives 
rise to £1 w 10~® while (72 is supposed to minimize 62 
with £2 ~ 10~^. It seems that, compared to three-level 
system, higher frequencies and amplitudes are needed to 
reach high fidelity of NOT-gate. In three-level system the 
iterative optimization algorithm is applied at most 5000 
times to reach such fidelities while with five levels 15000 
iterations were needed to obtain the results shown in fig- 
ures [TT] and [12] The leakage out of the qubit manifold 
would be the reason for this. 



VII. CONCLUSIONS 

In this paper we have shown that it is possible to opti- 
mize single-qubit gates for Josephson phase qubits by em- 
ploying quantum optimal control theory. We have con- 
sidered the realistic situation in which, in addition to the 
two computational basis states |0) and |1), higher energy 
states are present, which may lead to leakage. Typically 
microwave pulses with Gaussian modulation are used to 
induce transition between states |0) and yielding a 
quite high fidelity for long pulse durations. For the sake 
of definiteness. here we have focused on the NOT-gate 
single-qubit operation and searched for modulations of 
microwave pulses which optimize such operation, espe- 
cially for short-duration pulses. The numerical results 
obtained for a three-level system, and neglecting off- 
resonance terms, demonstrate up to ten orders of mag- 
nitude improvement in fidelity of a NOT-gate operation 
with respect to those obtained through Gaussian modu- 
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lations. To test the effect of possible imperfections in the 
pulses shape, we have studied the behavior of the fidelity 
as a function of the bandwidth of the pulse generator and 
showed that frequencies not bigger 2 GHz are needed to 
gain up to four orders of magnitude improvement. More- 
over, we have shown that the off- resonance elements of 
the Hamiltonian, which arc usually neglected, can be im- 
portant for optimized pulses, especially for short pulse 
duration times, due to the very high fidelity reached. 
We have also addressed the effect of the presence of a 
capacitively-coupled second qubit and showed that, even 
though the optimized pulses are obtained for a single 
qubit, they still lead to a high fidelity for a NOT-gate 
(up to two orders of magnitude improvement) especially 
for very short pulses. Finally, we were able to obtain op- 
timized pulses for a system with 5 energy levels in which 
the leakage outside qubit manifold is more severe. 

In conclusion, the two-interacting-qubit system de- 
serves for sure further attention. On the one hand, in 
order to improve the fidelity of a singlc-qubit operation, 
in the presence of capacitivc coupling, it seems that a 
way to switch the interaction on and off should be found 
even for optimized pulses of an isolate qubit. On the 
other hand, obtaining optimized pulses while including 
the interaction, would be a potential theoretical work to 
be done. 



1.0 

t (27i/S(o) 



FIG. 13: (Color on line) Examples of final optimized modu- 
lation of pulses (dashed lines), obtained from minimizing ei 
(top panel) and 62 (bottom panel), and the corresponding 
pulse with Gaussian modulation (solid lines) used as initial 
guess in optimization process with duration time T = 2 |^ in 
a system with five energy states. 
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